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Abstract 

We employ a novel new approach to study local quark-hadron duality using "truncated" moments, 
or integrals of structure functions over restricted regions of x, to determine the degree to which 
individual resonance regions are dominated by leading twist. Because truncated moments obey the 
same Q 2 evolution equations as the leading twist parton distributions, this approach makes possible 
for the first time a description of resonance region data and the phenomenon of quark-hadron duality 
directly from QCD. 
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I. INTRODUCTION 



The structure and interactions of hadrons at intermediate energies represents one of the 
most outstanding problems in the standard model of particle and nuclear physics. Many 
hadronic observables can be described at low energies in terms of effective, hadronic (me- 
son & baryon) degrees of freedom, while at high energies perturbative QCD has proved a 
highly successful approach to describing phenomena in terms of elementary quark and gluon 
constituents. The precise nature of the transition between the two regimes has remained 
shrouded in mystery, however, and represents a fundamental challenge to our understanding 
of strong nuclear interactions within QCD. 

One of the most intriguing connections between the low and high energy realms is the 
phenomenon of quark-hadron duality, in which one finds in certain cases dual descriptions 
of observables in terms of either explicit quark degrees of freedom, or as averages over 
hadronic variables. A spectacular example of such a duality is in inclusive electron-nucleon 
scattering. First observed by Bloom & Gilman in the early days of deep inelastic scattering 
(DIS) measurements [lj], this duality manifests itself in the similarity of structure functions 
averaged over the resonance region (which is characterized by hadronic bound states) and 
the scaling or leading twist function describing the high energy, deep inelastic continuum 
(characterized by scattering from free quarks) [2j. 

Unraveling the origin of the "Bloom-Gilman" duality from first principles has proved to 
be a major challenge in QCD. Until now the only rigorous connection with the fundamen- 
tal theory has been within the operator product or "twist" expansion, in which moments 
of structure functions are expanded as a series in inverse powers of the virtuality Q 2 of 
the exchanged photon. The leading, 0(1) term is given by matrix elements of (leading 
twist) quark-gluon bilocal operators, and is associated with free quark scattering, while the 
0(1/Q 2 ) and higher terms correspond to nonperturbative (higher twist) quark-gluon inter- 
actions. Bloom-Gilman duality is then interpreted in this language as the suppression of 
higher twist contributions to the moments Q|. 

Recent experimental data J| suggest that not only moments, but structure functions in 
individual resonance regions, such as the A, Sn or F 15 regions, closely resemble the leading 
twist structure functions over the same intervals. This indicates that duality also exists 
locally, in restricted regions of hadronic final state mass W. The appearance of this local 
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duality cannot, however, be explained with the theoretical tools presently at our disposal, 
and insight into the workings of duality for individual resonance regions has been confined 
to QCD-inspired models of the nucleon. As such our understanding of quark-hadron duality 
in nucleon structure functions within QCD is incomplete. 

In this paper we present a new approach to the study of local quark-hadron duality within 
a perturbative QCD context, using "truncated" moments of structure functions. The virtue 
of truncated moments is that they obey a similar set of Q 2 evolution equations as those for 

nn 

parton distributions [5|, [6|, which therefore enables a rigorous connection to be made between 
local duality and QCD. It allows us to quantify for the first time the higher twist content 
of various resonance regions, and determine the degree to which individual resonances are 
dominated by leading twist. 

Truncated moments were introduced several years ago by Forte et al. [s| to study struc- 
ture function moments for which small- £ data were not available. By restricting or truncating 
the integration region to some minimum value of the Bjorken x variable, one could avoid 
the problem of extrapolating parton distributions into unmeasured regions at small x. Later 
Kotlorz & Kotlorz [6|| developed an alternative formulation of the evolution equations which 
avoids the problem of mixing of higher truncated moments when evolving in Q 2 . 

In this work we partially follow the latter approach and apply it to the study of structure 
functions in the large-x region, populated by nucleon resonances. In particular, we study the 
Q 2 evolution of structure functions integrated over specific nucleon resonance regions. To 
facilitate such an analysis requires extension of the definition of the truncated moments to 
include both upper and lower truncations. We show that these "doubly truncated" moments 
also obey the same Q 2 evolution equations. Using recent high-precision data on the proton F 2 
structure function from Jefferson Lab and elsewhere, we quantify the size of the higher twists 
for the lowest three moments in various regions of W. This represents the first quantitative 
test of local duality in structure functions within a QCD framework. 

This paper is organized as follows. In Sec. [Til we review the essential elements of Q 2 
evolution via the DGLAP equations, and introduce truncated moments together with their 
evolution. We test the accuracy of our numerical evolution procedure in Sec. lIIIl and further 
study recent proton structure function data in the nucleon resonance region at W < 2 GeV. 
We divide the data into the three traditional resonance regions and extract the leading and 
higher twist content of each region. Finally, in Sec. [TV]we summarize the conclusions of this 
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analysis and outline future work. 



II. TRUNCATED MOMENTS AND EVOLUTION 
A. QCD Evolution Equations 

The Q 2 dependence of a parton distribution function (PDF) q(x, Q 2 ) is described in 
perturbative QCD (pQCD) by the DGLAP evolution equations [7||: 

dq(x,Q 2 ) a s (Q 2 ^ 



dt 27T 



'P <8> q) (x, Q 2 ) , (1) 



where t = In (Q 2 /Aq CD ), with Aqcd the QCD scale parameter, and the symbol <g> denotes 
the Mellin convolution, 

(P ® q) (*, Q 2 ) = f^p(-^ s (Q 2 )) q (y, Q 2 ) , (2) 

j x y \ y j 

between the parton distribution q and the splitting function (or the evolution kernel) P. 
In pQCD the latter can be expanded as a series in the strong running coupling constant 
a s{Q 2 )- For the nonsinglet (NS) case, q is one of the flavor nonsinglet combinations of 
quark distributions and P the corresponding NS splitting function. For the singlet case, on 
the other hand, q is a vector whose components are the flavor singlet combination of quark 
distributions and the gluon distribution, and correspondingly P is a 2 x 2 matrix of splitting 
functions. 

Taking moments, the convolution in Eq. (DD) turns into an ordinary product, and the 
evolution equations become ordinary first order differential equations in moment space n, 

dM n (Q 2 ) a s (Q 2 " 



dt 2tt 



ln(Q Z ) M n (Q ) , (3) 



which can be solved analytically. Here the n-th full moment of the parton distribution is 
defined as: 

M n (Q 2 ) = [ dxx 11 - 1 q(x,Q 2 ) , (4) 
Jo 

and the anomalous dimension, 

7n(Q 2 ) = f dzz n - 1 P{z,a s {Q 2 )) , (5) 
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is the moment of the splitting function P (z, as (Q 2 ))- The PDF can then be determined 
via the inverse Mellin transform, 

i rc+ioc 

Q&Q 2 ) = / dnx~ n M n (Q 2 ) . (6) 

From the definition in Eq. Q, the full moments are obtained by integrating the PDF 
over all values of the Bjorken variable, < x < 1. Since x is related to the invariant mass 
squared W 2 of the virtual photon-hadron system, W 2 = M 2 + Q 2 (1 — x) /x, where M is 
the nucleon mass, to reach the x — > limit requires infinite energy; hence in practice some 
extrapolation to x = is always needed to evaluate the moment. Similarly, at finite Q 2 one 
usually excludes from leading twist analyses the W < 2 GeV region in order to avoid \ow-W 
nucleon resonances, so an analogous extrapolation to x = 1 is also performed. 



B. Truncated Moments 



An alternative approach, which avoids uncertainties from unmeasured regions at low 

n 

and high x, makes use of the so-called "truncated moments" [5fl. In analogy with the full 
moments, the truncated moments of a PDF q(x, Q 2 ) are defined as: 



M n (x ,l,Q 2 ) = C dxx n ~ x q(x,Q 2 ) , 



■>■<) 



where the integration is restricted to x < x < 1 (the first two arguments in M. n denote the 
lower and upper limits of the integration). From the evolution equation jTJ, one can verify 
that the truncated moments satisfy: 

dM n (x ,l,Q 2 ) a s (Q 2 )_ rl 



dt 2tt 
where 



f dyy n ~ 1 q {y,Q 2 )G n ( X -^,Q 2 ) , (8) 



G n (x,Q 2 ) = f dzz n - 1 P{z,a s {Q 2 )) 

J X 



(9) 



is the truncated anomalous dimension. For x = 0, the latter reduces to the usual x- 
independent anomalous dimension, G n (0,Q 2 ) = 7 n (Q 2 ), which can be taken outside the 
integral in Eq. (jHJ). The right hand side then depends only on the n-th moment, and the full 
moments of PDFs evolve independently of each other. 

For nonzero xq, the residual y dependence in the truncated anomalous dimension leads 
to evolutions equations which are not diagonal in n. This can be seen by expanding 



G n (xo/y,Q 2 ) as a Taylor series around y — 1, and truncating the expansion at a finite 
order m. Accordingly, Eq. (jHJ) then turns into a system of coupled evolution equations: 



fc=0 



where 



<A*o) = Y, fc!(p _ fc)! and ^ = « J 



(11) 



y=X 



p=k 

Unlike the full moments, the evolution of the truncated moment of order n is determined 
by all truncated moments of order n + k, with k > 0. However, the series of couplings to 
higher moments converges, and can be truncated to any desired accuracy. One can solve 
Eq. (fTOl) to arbitrarily high accuracy by using a sufficiently large basis of truncated moments. 
For example, the higher moments {n > 2) can be calculated with excellent accuracy even for 
a small (m = 4) number of terms in the expansion of the truncated anomalous dimension. 
The first moment, on the other hand, is more sensitive to the truncation point Xq and the 
convergence of the truncated anomalous dimension for n — 1 is weaker than for the higher 
moments. 



C. Diagonal Formulation of Truncated Moments 

The evolution equations satisfied by the truncated moments can be formulated in an 
alternative way which avoids the problem of mixing of lower moments with higher moments 
0| . Inverting the order of integration on the right hand side of Eq. (jHJ) and introducing a 
new variable u = xo(y/x), the integral can be written as: 

/ dx i"" 1 (P <g> q) (x, Q 2 ) = (P; t ® M n ) {x , Q 2 ) , (12) 

Jx 

with 

P' n {z,a s {Q 2 )) = z n P(z,a s (Q 2 )) . (13) 

The evolution equation for the truncated moments then becomes: 

dM n (xo,l,Q 2 ) a s (Q 2 ) ,p, KA w n2 , , . 
j t = 2n ( p n ® M n ) {Xq, Q ) , (14) 
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which is very similar to the original evolution equation (CQ) for the PDFs. Here P' n plays the 
role of the splitting function for the truncated moments. The truncated moments therefore 
satisfy DGLAP evolution with a modified splitting function P(z, a>s{Q 2 )) — > z n P(z, as{Q 2 )) 
in the Mellin convolution. The advantage of this approach is that it can be successfully 
applied to any n-th moment and for every truncation point < x < 1, without the 
complication of mixing with higher moments. 

The evolution equations for the truncated moments can also be generalized to any subset 
in the x-region, x m \ n <x< x max . Writing the "doubly-truncated" moment of the PDF as: 

r 3^ max 

Mn(x m m,x mSLX ,Q 2 ) = / dx x n_1 q(x, Q 2 ) , (15) 

3- min 

its Q 2 evolution can be obtained by subtracting the solutions of truncated moments at the 
points x mhl and x max : 

where M n (x min , 1, Q 2 ) and A / l n (x max , 1, Q 2 ) both satisfy Eq. (PT4l) . 

III. DATA ANALYSIS 

The central aim of this study is to determine the extent to which nucleon structure 
function data in specific regions in x (or W) are dominated by leading twist. This can 
be done by constructing empirical truncated moments and evolving them to a different 
Q 2 using one of two methods. Namely, (i) the structure functions are evolved and the 
corresponding truncated moments the calculated, or (ii) the moments are evolved directly 
using the evolution equations in Eq. (TT4l) above. We found the results of both methods to be 
essentially equivalent. In the study of the proton structure function data there is, however, 
difficulty in applying the target mass corrections (TMCs) using the latter method. Here 
in principle one can derive and solve the evolution equations for the target mass corrected 
moments (the so-called truncated Nachtmann moments), which contain the TMCs explicitly. 
In practice, to avoid this problem we shall utilize method (i): we evolve the structure 
functions, correct them for TMCs, and finally calculate their moments. 

Deviations of the evolved moments, computed to next-to-leading order (NLO) accuracy, 
from the experimental data at the new Q 2 then reveal any higher twist contributions in the 
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original data. In particular, we will analyze recent data on the proton F% structure function 
from Jefferson Lab covering a range in Q 2 from ^ 1 GeV 2 to ~ 6 GeV 2 . 

The evolution of the measured truncated moments requires the structure function to be 
decomposed into its nonsinglet and singlet components. Without performing a global pQCD 
analysis of the structure function data, it is a priori unknown which parts of the structure 
function are singlet and nonsinglet. To proceed, we shall assume that in our region of 
interest, at moderate to large x, the proton structure function is well approximated by its 
nonsinglet component, and will evolve the truncated moments as nonsinglets. The accuracy 
of this approximation will improve with increasing order of the truncated moments. 



A. Evolution of Truncated Moments 



We test the accuracy of the nonsinglet evolution by first evolving a trial structure function 
whose decomposition into its nonsinglet and singlet components is known. The trial function 
is evolved exactly, with its nonsinglet and singlet components computed separately, and 
also evolved under the assumption that the total function can be treated as a nonsinglet. 
A comparison of the discrepancy between the two evolved truncated moments can then 
reveal the accuracy of the nonsinglet evolution of the various moments as a function of the 
truncation region. 

There are many methods to solve the Q 2 evolution equations for the truncated moments 
8|. The simplest and most direct is to solve the equations by brute force using a suitable 
numerical integration routine. In this work we use the method of Ref. [s| for the evolution. 

To illustrate the method of direct moment evolution, we consider the evolution of the non- 
singlet truncated moment, A4^ s , to leading order in as (although in practice our numerical 
results are performed at NLO): 



l dy 



x 



-» _ p 



dr J x y \y 

where the leading order NS splitting function is 



,(0) 



NS 



.r 



- MT(y,i,r) 



l + z 2 



:S(1 



and instead of t we use the variable r, where 



as(Q 2 
<*s{Ql)\ 



(17) 



(18) 



(19) 




Figure 1: Ratio of the truncated moments of F^ evolved from Q 2 = 25 to 1 GeV 2 , using NS and 
full evolution, versus the truncation point H^ max (or x m - m ), for the n = 2 (solid), 4 (dotted) and 6 
(dashed) moments. 

with cts(Q 2 ) = 4:7i / [p Q ln(Q 2 /AQ CD )] the running coupling constant at leading order, and 
Po = 11 — 2Nf/3 for Nf quark flavors. By dividing the variables r and x into small steps, 
the evolution from Tj to Tj+\ can be written as: 

+ ^E^(5)"^(5)^fe,l,r 3 ), (20) 

k=i N v ' 

where At, = r J+1 — r, and Ax*. = — Xk-\ are the steps at positions j and k, and N x is 
the number of steps in x. The final truncated moment at tn t is then obtained by repeating 
the step in Eq. (N T - 1) times. 

For the trial function we take the leading twist proton F 2 structure function computed 
from the MRST2004 PDF fit flo| . The n = 2, 4 and 6 truncated moments of are then 
evolved from Q 2 = 25 GeV 2 to 1 GeV 2 using NS evolution, and compared with the exact 
results using singlet and nonsinglet evolution. The ratios of these are plotted in Fig. Q] as a 
function of the truncation point W max , where W^ ax = M 2 + Q 2 (l/x min — 1). 
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Generally the differences between the full and NS evolution are of the order 2-4% for 
1.2 < W m ax ~ 2 GeV, the traditional nucleon resonance region, and increase with increasing 
Wmax- Note that at Q 2 = 1 GeV 2 , W max = 2 GeV corresponds to z min = 0.24. For the n = 2 
moment, which is most sensitive to singlet evolution, the differences do not exceed ~ 4% 
for 1.2 < W max ^ 2 GeV. As expected, for the higher moments the differences are smaller, 
^ 2 — 3% for n = 4 and n = 6 for 1.3 < W max < 2 GeV. In the region relevant for our study 
one can therefore safely conclude that the error introduced by evolving the F\ moments 
as nonsinglets is less than 4%. This uncertainty will be included in the errors in our final 
results. 



F% data 



fit 



B. Extraction of Higher Twists 

Having tested the accuracy of the nonsinglet evolution we now turn our attention to 
the analysis of the F% structure function data. In Figj2| (top panel) we compare the F% 
resonance data from Jefferson Lab experiment E91-110 [ill] at Q 2 = 1 GeV 2 (triangles) with 
an empirical fit [2l| to the data (dashed), and with leading twist fits to the deep-inelastic 
20fl with (solid) and without (dot-dashed) target mass corrections. The resonance 
2lJ describes the Ff data to better than 3% over the range < Q 2 < 8 GeV 2 and W 2 
from the inelastic threshold up to 10 GeV 2 . 

Since the data at low Q 2 contain significant contributions arising from kinematical M 2 /Q 2 
corrections (which, although subleading in Q 2 , contribute at leading twist), a direct com- 
parison of data with the leading twist structure function requires the inclusion of TMCs. 



We d 
Ref. 



o so here by applying the standard TMC prescription for F 2 from Ref. 



121 ] (see also 



131 ] for a review of TMCs). As is evident from Fig. [2j the leading twist fit to the DIS 
data, including TMCs, agrees well with the average F| data in each resonance region. A 
comparison of this fit with the DIS data at Q 2 = 25 GeV 2 (bottom panel of Fig. [2]) shows 
the excellent agreement between the leading twist function and the data at this scale. 

More specifically, the comparison of the data with the target mass corrected leading 
twist function illustrates the intriguing phenomenon of Bloom-Gilman duality, where the 
data in the resonance region oscillate around, and on average are approximately equal to, 
the leading twist function This duality reveals itself in the relatively small value of the 
higher twist contributions at these scales observed in recent high-precision F| measurements 
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Figure 2: (Top panel) Comparison of proton F% data from JLab experiment E91-110 



1 GeV 2 (triangles), and the fit 



HI] at Q 2 



2l| to the data ("Resonance Fit", dashed), with a global fit of 



2fJ| . Shown is the leading twist DIS fit with 



DIS data used to extract target mass contributions 
("LT+TMC", solid) and without ("LT", dot-dashed) the TMCs. The vertical lines indicate the 
extent of the second (<Sii) resonance region. (Bottom panel) DIS data at Q 2 = 25 GeV 2 (circles) 
compared with the LT and LT+TMC fits. 



f. Q 



— see Ref. |2j for a review of the data. Note that the nonzero value of F 2 with TMCs in 
the limit x — > 1, which is related to the so-called threshold problem, introduces a small 
additional uncertainty into structure function analyses at low Q 2 fl^. 

To determine the extent to which the data at low Q 2 are dominated by leading twist, we 
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assume that the data beyond some large Q 2 value are dominated by twist-2 contributions. 
In view of the comparison with the data in Fig. El in this analysis we take this scale to 
be Q 2 = Ql = 25 GeV 2 . This assumption is also consistent with most global analyses 
's, which fit leading twist PDFs to structure function data down to Q 2 ~ 1-2 GeV 2 



of PD 
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17l |. Although these analyses typically exclude low-PF resonance data, in practice 
there is little contribution to the low moments from the resonance region W < 2 GeV. 
The analysis method then proceeds in four main steps: 

1. For each W 2 region of interest, the x range to be covered (Ax) is calculated at the 
particular (lower) Q 2 where the leading twist contribution is to be extracted. For the 
second (<Sii) resonance region, for example, this is indicated by the vertical lines in 
Fig. [2] (top panel). 



2. The structure function extracted from a precision fit to data [20|] at the starting scale 
Q 2 = 25 GeV 2 is evolved in leading twist down to lower Q 2 and the TMCs are applied. 
Note that at the higher Ql value the same interval Ax is used, as e.g. indicated by 
the vertical lines in Fig. [2] (bottom panel), which corresponds to a higher W range. 

3. The truncated moment of the evolved leading twist, target mass corrected structure 
function from Step 2 is calculated for the interval Ax defined in Step 1. 

4. The truncated moment of the resonance data at Q 2 is calculated in the interval Ax 
and compared with the result of Step 3. 

After evolving down to Q 2 = 1 GeV 2 and applying the target mass corrections, the n = 2 
truncated moment M. 2 is shown in Fig.[3](a) as a function of where it is compared with 

the moment of the actual data at Q 2 = 1 GeV 2 . The difference between the evolved curve 
and the data attests to the presence of higher twist contributions in the data at Q 2 = 1 GeV 2 . 
The importance of the TMCs is also clearly evident, and these in fact reduce the difference 
between the leading twist moment and the data by some 40% for large PF max . The ratio of 
the truncated moments of the data to the leading twist in Fig. [3](b) illustrates that without 
TMCs the leading twist moment differs from the data by ~ 20% for W mauX > 1.5 GeV. After 
correcting for TMCs, the size of the apparent higher twists is reduced to ~ 15%. It is 
imperative, therefore, that the kinematical effects associated with finite values ofQ 2 /v 2 be 
properly accounted for before drawing any conclusions about higher twists from data. 
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Figure 3: (a) Truncated moment M.2 as a function of the truncation point W max at Q 2 = 1 GeV 2 , 
evolved as leading twist (LT) from Qq = 25 GeV 2 with (solid) and without (dashed) target mass 
corrections (TMC), and compared with the moment calculated from data (dotted) at Q 2 = 1 GeV 2 . 
(b) Ratio of the M.2 truncated moments of the data to the leading twist + TMC (solid), and data 
to leading twist without TMC (dashed) at Q 2 = 1 GeV 2 . 

Note that the truncated moments displayed in Fig. [3] are computed over the range W t ^ < 
W < W max , where the W t h = M + m n is the inelastic threshold. This is consistent with the 
assumption that the truncated moments at Ql = 25 GeV 2 are entirely of twist-two, since 
the elastic cross section contributes only to the higher twist part of the structure function. 
For a meaningful comparison, we therefore do not include the elastic contribution at lower 
Q 2 . At Q 2 = 1 GeV 2 this corresponds to the integration range x m \ n < x < x t h, where 
x th = [1 + m,(m 7r + 2M)/Q 2 Y l ~ 0.78. 

The results in Fig. [3] give a clear indication of the magnitude and sign of higher 
twists in the data at Q 2 = 1 GeV 2 . To quantify the higher twist content of the spe- 
cific resonance regions, and at different values of Q 2 , we consider several intervals in 
W: W 2 h < W 2 < 1.9 GeV 2 , corresponding to the traditional A(1232) (or first) reso- 
nance region; 1.9 < W 2 < 2.5 GeV 2 for the Sii(1535) (or second) resonance region; and 
2.5 < W 2 < 3.1 GeV 2 for the _Fi5(1680) (or third) resonance region. The n = 2 truncated 
moments corresponding to these regions are plotted in Fig. [4] for various Q 2 values, from 
Q 2 = 1 GeV 2 to Q 2 = 6 GeV 2 . Shown are ratios of moments calculated from the data to the 
moments obtained from NLO evolution of the leading twist moments from Ql = 25 GeV 2 , 
corrected for target mass effects. Below Q 2 = 1 GeV 2 the applicability of a pQCD analysis 
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Figure 4: Q 2 dependence of the ratio of n = 2 truncated moments M.2 calculated from the data 
and from leading twist evolution from Q\ = 25 GeV 2 (including TMCs), for various intervals in 
W: the first (A) resonance region (diamonds), second (Sn) resonance region (squares), the first 
and second combined, corresponding to W^ax = 2.5 GeV 2 (dotted curve), third (-F15) resonance 
region (open circles), first three regions combined, W^ax = 3.1 GeV 2 (dashed curve), and the entire 
resonance region Wmax = 4 GeV 2 (filled circles). Note that some of the points are offset slightly for 
clarity. 

becomes doubtful and the decomposition into leading and higher twists is no longer reliable. 

The results indicate deviations from leading twist behavior of the entire resonance region 
data (filled circles in Fig. [4j at the level of ^ 15% for all values of Q 2 considered, with 
significant Q 2 dependence for Q 2 < 4 GeV 2 . This is made more explicit in Fig. El where the 
higher twist contributions to M.2 (defined as the difference between the total and leading 
twist moments) are shown as ratios of the moments evaluated from the data. 

The strong Q 2 dependence of the higher twists is evident here in the change of sign around 
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Figure 5: Q 2 dependence of the fractional higher twist (HT) contribution to the n = 2 truncated 
moment data, for various intervals in W (as in Fig. HJ. 

Q 2 = 2 GeV 2 , with the higher twists going from « -10% at Q 2 = 1 GeV 2 to « 10-15% for 
Q 2 « 5 GeV 2 . The slope at Q 2 « 1 - 2 GeV 2 would be decreased if the full NS + singlet 
evolution were performed, as evident from Fig. CD, since the NS-only evolution leads to a few 
percent overestimate of the LT+TMC results. At larger Q 2 the higher twists are naturally 
expected to decrease, once the leading twist component of the moments begins to dominate. 
Note that the extraction of higher twists beyond Q 2 = 6 GeV 2 would require evolution from 
a starting scale larger than the Ql = 25 GeV 2 used in this analysis. At larger Q 2 , however, 
data in the large-x region, which determines the behavior of the resonances after evolution to 
lower Q 2 , are not well determined, making extraction of higher twists beyond Q 2 ~ 6 GeV 2 
problematic at present. 

Turning to the individual resonance regions, the results in Figs. [Hand [5] show that in the 
A region (diamonds) the higher twist contributions are smallest in magnitude at large Q 2 , 
decreasing from ps —15% of the data at Q 2 = 1 GeV 2 to values consistent with zero (within 
errors) at larger Q 2 . The higher twists are largest, on the other hand, for the Sn region 
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(squares), where they vary between « -15% of the data at Q 2 = 1 GeV 2 and 20-25% at 
Q 2 ~ 5 GeV 2 . Combined, the higher twist contribution from the first two resonance regions 
(dotted curve) is ^ 15% in magnitude for all Q 2 . The rather dramatic difference between 
the A and the Su, may, at least in part, be due to the choice of the differentiation point of 
W 2 = 1.9 GeV 2 . A lower W 2 choice, for instance, would lower the higher twist content of 
the Sn at large Q 2 , while raising that of the A. However, our W 2 choice corresponds to the 
local minimum between these two resonances in the inclusive spectra, and is the one most 
widely utilized. 

The higher twist content of the Fi 5 region (open circles) is similar to the Su at low Q 2 , 
but decreases more rapidly for Q 2 > 3 GeV 2 . The higher twist content of the first three 
resonance regions combined (dashed curve) is ^ 15-20% in magnitude for Q 2 < 6 GeV 2 . 
Integrating up to = 4 GeV 2 (filled circles), the data on the n = 2 truncated moment 

are found to be leading twist dominated at the level of 85-90% over the entire Q 2 range. 

The results in Figs. 0] and [5] contain two sources of uncertainty: from the experimental 
uncertainty on the F 2 data (statistical and systematic) , and from the nonsinglet evolution of 
the data. For the experimental error we take an overall uncertainty of 2% for all truncated 
moment data, with the exception of the n = 4 and n = 6 moments for W 2 iax = 1.9 and 
4 GeV 2 , where the experimental uncertainties are 4% and 3% for A4 4 , and 5% and 4.5% 
for M.q for the two W regions, respectively. The evolution error is estimated by comparing 
the nonsinglet evolution with the full evolution using the MRST fit, as in Fig. [U with the 
appropriate correction factor applied at each Q 2 and W interval. We do not assign an 
uncertainty for the structure function at the Qq = 25 GeV 2 input scale, as this is negligible 
for all but the A region analysis at Q 2 = 6 GeV 2 , where we estimate it to be < 3%. 

Note that the size of the error bars on the Su and F 15 data points at the lowest Q 2 are 
larger than those for the A. This is due to the fact that these resonances lie at higher W, 
and hence at lower values of x compared to the A, and where the NS approximation to the 
full evolution is not as good as at large x. The error bar on the A is smaller at the same 
low Q 2 values since it appears at larger x, where the NS versus full evolution differences are 
smaller. 

The relatively small size of the higher twists at scales ~ few GeV 2 is consistent with 
the qualitative observations made in earlier data analyses about the approximate validity of 
Bloom-Gilman duality Q|. In this analysis we are able to for the first time quantify precisely 
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the degree to which this duality holds as a function of Q 2 (see also Ref. [18]). The fact that 



duality works better (higher twists are smal 
borne out in various quark model studies j3, 



erj when more resonances are included is also 
ISO. 
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Similar behavior is found also for the n = 4 and n = 6 truncated moment ratios, illus- 
trated in Figs. [6] and [3, respectively. For the higher moments, the overall magnitude of the 
higher twists is qualitatively similar to the n = 2 moments, although the Q 2 values at which 
they start decreasing in importance are larger. At low Q 2 values the higher twist contri- 
butions are also relatively larger for higher moments: at Q 2 = 1 GeV 2 , for example, the 
magnitude of the higher twist component of the W 2 < 4 GeV 2 region increases from ~ 10% 
for the n = 2 moment, to ~ 15-20% for n = 4, and ~ 20-30% for n — 6. This behavior 
can be understood from the relatively greater role played by the nucleon resonances and the 
large- x region, which is emphasized more by the (x-weighted) higher moments. 

IV. CONCLUSION 

Quark-hadron duality in nucleon structure functions remains an intriguing empirical phe- 
nomenon which challenges our understanding of strong interaction dynamics, as one transi- 
tions from hadronic degrees of freedom in the nucleon resonance region to quarks and gluons 
in the deep inelastic continuum. Until now the only rigorous connection with QCD has been 
for moments of structure functions analyzed within the twist expansion. Any insight about 
the workings of duality for individual resonances, or specific resonance regions, has been 
confined to QCD-inspired models of the nucleon. 

In this paper we have presented a new approach to the study of local quark-hadron 
duality within a perturbative QCD context, using so-called truncated moments of structure 
functions. The fact that truncated moments obey a similar set of Q 2 evolution equations 
to the DGLAP equations for parton distributions enables a rigorous connection to be made 
between local quark-hadron duality and QCD. It allows us to quantify for the first time 
the higher twist content of various resonance regions, and determine the degree to which 
individual resonance regions are dominated by leading twist. 

We find deviations from leading twist behavior of the truncated moments of the resonance 
region data (W < 2 GeV) at the level of < 15% for Q 2 > 1 GeV 2 . Significant Q 2 dependence 
in the ratio of moments of data to leading twist is evident for Q 2 < 3 GeV 2 , with the higher 
twists changing sign around Q 2 = 2 GeV 2 . For the n = 2 truncated moment, M.2, the higher 
twists are found to vary from « -10% at Q 2 = 1 GeV 2 to « 10-15% at Q 2 » 5 GeV 2 . 

Separating the W < 2 GeV data into the three traditional resonance regions, our results 
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indicate that at a scale of Q 2 = 1 GeV 2 the A resonance region contains about —15% higher 
twist component of the total M. 2 , but is consistent with zero at larger Q 2 . The higher twists 
in the second (Sn) and third (-F15) resonance regions are larger in magnitude, with the Sn 
ranging from « -15% at Q 2 = 1 GeV 2 to 20-25% at Q 2 ~ 5 GeV 2 , and the F 15 varying 
from and 15% over the same range. 

Similar behavior is found also for the n = 4 and n = 6 truncated moments. Here 
the relatively greater role played by the resonances due to the large-x enhancement leads 
to larger higher twists at the same Q 2 . At Q 2 = 1 GeV 2 , for example, the higher twist 
component of the W < 2 GeV region increases from around —10% for M. 2 to ~ —15% for 
Mi, and ~ -25% for Me- 

In contrast to earlier analyses of duality using complete moments of structure functions 
which have quantified the total higher twist content over all x, this analysis in terms of 
truncated moments reveals the distribution of higher twist corrections over various regions 
in x (or W). Note that, unlike many previous moment analyses, an effort was also made to 
quantify the uncertainty associated with evolving the structure function data as a nonsinglet, 
which was found to be ^ 4%. 

While this analysis has been to some extent exploratory, it has illustrated an encourag- 
ing new approach to quantifying and understanding local Bloom-Gilman duality within a 
well-defined theoretical framework. It opens the way to further study of local duality in 
other structure functions, such as the longitudinal structure function Fi or spin-dependent 
structure functions. 
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